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SUMMARY 


A  study  of  the  strength  of  materials  under  large  deformation  has 
been  conducted.  The  sample  to  be  tested  was  in  the  form  of  a  thin  disk  and 
was  crushed  between  a  fixed  and  a  moving  anvil.  Velocities  ranging  from 
10“3  to  10^  cm/sec  were  obtained  by  driving  the  anvil  either  by  a  testing 
machine  or  by  a  compressed-air  gun. 

A  simplified  analysis  of  the  problem  is  given,  and  the  assumptions 
used  are  discussed  and  compared  with  the  experimental  data.  The  relation¬ 
ship  between  internal  energy  of  the  material  and  deformation  is  found  to  be 
U  =  'f  In  ( Zq/z )  where  U  is  energy  per  unit  mass,  zQ  and  z  are  initial  and 
final  thickness  of  the  disk  and  T  is  a  constant,  the  "strength"  of  the 
material.  Variations  in  T  are  discussed  in  terms  of  strain-rate  effects  and 
strain-hardening  effects. 

The  materials  tested  were  copper,  99.99  per  cent  pure  aluminum, 
6061-T4  and  -T6  aluminum.  Nylon.  Lexan,  Teflon,  polyethylene  and  poly¬ 
propylene. 
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1.  INTRODUCTION 


The  purpose  of  research  dene  under  this  contract  is  to  develop 
a  theoretical  and  experimental  basis  for  predicting  the  behavior  of  materials 
under  Impact.  The  materials  of  interest  range  from  simple,  homogeneous 
materials  such  as  metals  or  ceramics  to  complex  structures  made  up  of  a 
variety  of  materials.  Impact  velocities  of  interest  range  from  a  f~w  hundred 
to  a  few  thousand  feet  per  second. 

One  approach  to  solving  the  problem  has  been  Lo  develop  a  versatile 
computer  program  to  solve  the  transient  flow  field  for  a  general  axisymmetric 
impact  problem.  Past  work  related  to  this^as  resulted  in  computer  cedes 
suitable  for  compressible  inviscid  fluids.  A  Work  is  in  progress  at  several 
centers  tc  include  viscosity  of  the  material  in  the  program,  but  so  far,  only 
an  artificial  viscosity  applicable  to  a  Newtonian  fluid  has  been  used,  and 
success  is  uncertain.  ^ 

In  the  velocity  range  of  interest  here,  strength  and  plastic-flow 
effects  are  as  important  as  inertial  effects  and  probably  dominate.  For  this 
reason,  our  approach  has  used  elasticity  and  plasticity  theory  as  well  as  a 
fluid-dynamics  approach.  Tensor  formulation  of  finite-deformation  problems 
has  been  explored  to  attempt  to  find  the  proper  mathematical  framework  for  the 
problem.  A  preliminary  computer  code  has  been  developed  for  a  simplified 
problem  consisting  of  the  edge-on  impact  of  two  finite-thickness  sheets 
having  infinite  length  and  width.  This  work  has  been  reported.  ^ 

A  second  approach  to  the  general  problem  has  been  to  formulate 
simple  models  of  the  impact  process  based  on  measurable  dynamic  material 
properties  then  to  develop  correlation  formulae  relating  observed  results  to 
measured  properties. 


1R.  L.  Bjork,  Effects  of  a  Meteoroid  Impact  on  Steel  and  Aluminum 
in  Space,  Technical  Report  P-1662,  Rand  Corp.  ,  Santa  Monica,  Calif.,  Dec. 
16,  1958. 
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T.  D.  Riney,  "Visco- Plastic  Solution  of  Hypervelocity  Impact 
Cratering  Phenomena, "  Proc.  of  Sixth  Symposium  on  Hypervelocity  Impact, 
Vol.  II,  p.  105,  Firestone  Rubber  Co,,  Contract  DA  31-124-ARO(D)-16, 
Cleveland,  Ohio,  Aug.  1963. 

^Technical  Reports  on  Contract  DA  19-129-AMC-150(X),  Final 
Report  on  Phases  I  and  II,  28  June  1964,  Semiannual  Report,  Phases  III,  IV, 
and  V,  November  1964, 
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In  this  report  period,  tne  work  has  been  aimed  at  devisinq  and 
analyzing  a  suitable  experiment  for  measuring  manorial  strength  under  large 
deformation  or  flow.  This  information  is  essential  for  progress  to  be  made  in 
theory  or  for  empirical  comparisons  of  materials.  The  experiment  chosen  is 
that  of  crushing  a  disk  between  two  moving  rigid  anvils.  Reasonable  simplify¬ 
ing  assumptions  make  an  analysis  possible,  and  the  necessary  experimental 
measurements  can  be  made.  The  results  of  this  work  are  the  subject  of  this 
report. 
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2.  DISK-CRUSHING  EXPERIMENT 


The  disk-crushing  experiment  uses  two  hardened-steel  anvils,  one 
fixed  and  one  driven  by  a  testing  machine  or  shot  from  a  gun.  The  disk  is 
placed  on  the  fixed  anvil  and  some  lubricant  used  to  reduce  friction.  If  a 
relatively  thin  disk  is  crushed  between  two  unyielding  anvils  at  a  velocity 
which  is  low  compared  with  wave  velocity  in  the  material,  we  may  assume  that 
the  material  is  incompressible,  inelastic,  und  that  flow  is  laminar  without 
friction  with  the  anvils.  With  unyielding  anvils  and  no  stored  elastic  energy, 
the  energy  of  crushing  can  be  determined  by  measuring  the  testing  machine 
force  and  travel  distance  or  by  measuring  velocity  and  mass  of  the  anvil  shot 
from  a  gun.  The  time  history  of  the  deformation  can  be  determined  as  well  as 
the  total  energy  involved.  In  practice,  the  assumptions  prove  to  be  justified 
or  corrections  can  be  made  over  limited  ranges  of  the  variables.  The  analysis 
of  the  problem  and  experimental  procedures  and  results  will  be  given  in  the 
following  sections. 


2 . 1  Differential  Equations  of  Disk-Crushing 
The  geometry  of  the  problem  is  as  shown. 
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The  subscripts  used  are  as  follows:  o  =  initial  condition,  f  =  final  condition, 
s  =  measurement  made  at  the  moving  surface. 

Neglecting  body  forces,  the  equations  of  continuity,  motion,  and 
energy  for  the  system  may  be  written  in  cylindrical  coordinates  as  follows: 
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Motion  (r  component) 
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The  7"  terms  are  the  components  of  the  stress  tensor. 
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These  equations  may  be  simplified  by  the  symmetry  of  the  problem  and  by  the 
assumptions  mentioned. 

One  simple  flow  pattern  assumes  that  the  z  velocity  decreases 
linearly  from  that  at  the  moving  anvil  to  zero  at  the  fixed  anvil.  We  assume 
3p/<9z  =  O  during  the  time  of  interest  after  the  first  waves  have  passed 
through  the  disk.  If  the  material  is  incompressible,  the  radial  flow  is  then 
determined.  The  rate  at  which  the  piston  displaces  material  at  a  portion  of 
the  circular  disk  of  radius  r  is  equal  to  the  rate  at  which  material  flows  out 
through  the  curved  sides  of  the  cylinder  at  r. 


Vr3ui  =  arr  (z0  -  zs~)  wr 

r  ir 

y; ;  - - - - 

*  (Zo-^s) 


(6) 


Vs  (z.-z) 

(Zo  ~zs) 


z 


Zs  (7) 


vs  is  a  function  of  time  and  describes  the  motion  of  the  moving  anvil. 

By  symmetry  =  0  ,  V~e  =  0  ,  Tr0  -  0.7'zg  =  O 

The  continuity  equation  (1),  becomes 
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Equations  (6)  and  (7)  satisfy  this  continuity  equation. 

The  0  component  equation  of  motion  is  identically  zero. 
The  r  and  z  component  equations  become: 


P  f  i r^t)  --  $£_{  s=^r  +  72r_  _  Zk.  /  o) 

'  K  dt  r  Br /  dr  [  dr  r  r  I 


I 


'jlTz  ,  - 

at  2 


Tzz 

c?Z 


(10) 


,0 


Considering  these  equations  with  (6)  and  (7),  we  see  that  by 
observing  the  dynamic  deformation,  vs  and  zQ-  zs,  we  can  make  important 
deductions  as  to  material  properties.  By  the  assumed  flow,  relative  material 
motion  in  the  z  direction  is  uniform  throughout  the  disk,  so  /<£z  —  Q 

vz  and  vr  can  then  be  determined  by  observing  vs  and  zs.  Some  prelimina  y 
wortc  in  analyzing  the  equations  of  motion  is  reported  in  the  experimental 
section. 


In  the  energy  equation,  the  pressure  term  is  zero,  as  is  seen  by 
substituting  Equations  (6)  and  (7)  into  it.  This  is  a  result  of  ^he  assumption 
of  incompressibility.  Heat  transfer  is  neglected  so  the  q  terms  are  zero. 

The  energy  equation  then  becomes 
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This  equation  can  be  simplified  further  by  using  Equations  (6)  and 
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Note  that  "V  is  opposite  in  sign  to  7V  and  and 

c^vz/cPz  Is  opposite  tooivj/c^r  and  vxfi.  The  term  in  parenthesis  may  be 
lumped  into  one  strength  term  7V 


Equation  (12)  provides  a  means  of  measuring  T~ and  also  ^T" n  / 
T qq  and  'T'zz  since  when  the  material  is  deforming  to  a  large  degree,  the 
three  terms  are  probably  about  equal. 


At  low  impact  velocities,  the  strain  is  dominated  by  the  dissipative 
term  and  is  probably  fairly  uniform  throughout  the  material.  Under  the 
assumptions,  the  dissipation  is  uniform  throughout  the  disk  so  the  substantial 
derivative  can  be  replaced  by  dU/dt.  Thus,  the  energy  equation  is  of  the 
simple  form 


-6- 


dU 

dT 
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where  zQ  and  Zf  are  initial  and  final  thickness  of  the  disk. 

The  actual  nature  of  'T' can  be  investigated  by  measuring  deviations 
from  Equation  (13)  as  strain  rate  or  total  strain  is  changed.  Various  useful 
approximate  equations  can  be  written  based  on  experimental  evidence.  As 
an  example,  7“  may  be  written  with  terms  giving  the  strain-rate  effect,  the 
work-hardening  effect  and  possible  work  softening  effects. 

-  z/z  -  U  _  u 
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The  a  coefficient  is  determined  b}  static  strength,  the  b  term  is  determined 
by  the  magnitude  of  the  strain-rate  effect  while  c  determines  the  strain-rate 
at  which  the  effect  becomes  significant.  Similarly,  d  and  f  determine  the 
work-hardening  effect  depending  on  total  energy  dissipated  or  U  and  g  and  h 
represent  a  work- softening  effect  which  dominates  at  higher  energies. 

The  functional  relationship  among  these  variables  may  be  quite 
different  from  that  indicated  here  for  s^me  materials  and  some  impact  conditions. 
Under  the  brief  exploratory  program  conducted  so  far,  using  copper,  aluminum, 
and  some  plastics,  the  ideas  embodied  in  Equation  (14)  seem  adequate  since 
7~  is  only  weakly  dependent  on  strain  or  strain  rate  and  the  constant  term,  a, 
dominates. 


The  experimental  difficulties  of  separating  friction  effects  from  strain- 
hardening  effects  make  the  determination  of  d  and  f  difficult. 

One  approach  to  correcting  for  friction  uses  the  following  model. 

An  infinitesimal  area  2TTx  dr  on  the  disk  at  radius  r  from  the  center  moves  with 
velocity  r  vs/2(z0-zs).  The  power  going  into  heat  due  to  I. iction  on  this  area 
is 

h  Cf  7T  \r5  cLr 
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where  F  is  the  force  per  unit  area  under  the  test  conditions  and  Cj  is  the 
coefficient  of  friction.  Integrating  from  r  =  O  to  the  outer  radius  R  and 
multiplying  by  2  to  include  both  faces  of  the  disk;  gives  the  total  power 
going  into  friction  as 

p  Cf  7T  Vs  R  3 _ 

f  3  (Z„-  2=^ 

vs/{zQ-zs)  can  be  expiessed  in  terms  of  R  using  Equation  (6\*  Pf  then  becomes 

Pf  -  F  Cf  T  LTr  R- 

or 
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where  W  is  the  frictional  energy.  Integrating  from  initial  to  final  disk  radius 

3 

Disk  mass  is  proportional  tc  R  for  disks  of  the  same  relative 
dimensions,  so  energy  per  unit  mass  contributed  by  friction  is  independent  of 
disk  size  for  disks  of  the  same  relative  length  to  diameter. 

Experiments  have  been  conducted  in  a  testing  machine  in  which  a 
disk  is  compressed  in  stages  and  comparisons  made  with  new  d,sks  (not  work 
hardened)  corresponding  to  each  of  the  stages.  These  tests  clearlv  snow  the 
separation  of  effects  of  work  hardening  and  friction  but  have  not  yet  been 
analyzed  to  obtain  quantitative  values  for  the  coefficients  in  Equations  (14) 
and  (IS). 


2  „  2  Experiments i.  Work 

Disks  made  of  copper,  aluminum, and  various  plastics  have  been 
deformed  by  two  methods:  (1)  by  impacting  them  at  velocities  of  ten  to  seventy 
meters  p^r  second  with  a  hardened- steel  cylinder,  and  (2)  by  compressing  them 
in  a  testing  machine  under  static  conditions.  The  deformation  ra^e  in  experi¬ 
ments  using  these  two  methods  varied  by  as  much  as  a  factor  of  10^.  It  was 
anticipated  that  such  a  variation  would  show  any  strain-rate  dependency  in  the 
materials  tested  merely  by  observing  variations  in  using  Equation  (13). 
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Most  of  the  uisks  were  of  the  same  original  size;  1/4  inch  diameter 
by  1/8  inch  thick..  However,  copper  disks  of  other  sizes  were  also  used. 
Materials  used  were  anrealed  and  unannealed  copper,  99,99  per  cent  pure 
aluminum,  6061-T4  and  b061-T6  aluminum  alloy.  Nylon,  Lexan,  Teflon, 
polyethylene  and  polypropylene. 

The  specific  energy  (energy  per  gram  of  material)  and  the  degree  of 
deformation  (the  ratio  of  the  final  disk  thickness  to  its  original  thickness) 
have  been  determined  for  each  impact  shot  and  for  each  compression  test. 

In  Figures  1  to  7 ,  deformation  is  plotted  on  a  logarithmic  scale  versus  specific 
energy.  Figure  1  shews  data  for  copper,  'igure  2  shows  data  for  aluminum, 
and  Figures  3  to  7  show  data  for  plastics.  If  the  relationship  expressed  in 
Equation  (13)  is  correct,  lines  drawn  through  the  points  in  these  figures  will 
be  straight,  indicating  that  specific  energy  is  directly  proportional  to  the 
logarithm  of  the  final  and  original  disk  thickness  ratio.  This  is  seen  to  be 
approximately  true  for  all  cases. 

In  Figure  1,  the  data  for  copper  show  that  there  is  no  observable 
difference  between  disks  of  different  sizes.  This  indicates  that  friction 
losses  are  proportional  to  mass  for  disks  of  the  same  relative  dimensions. 

The  best- fit  curves  to  the  copper  data  are  shown  in  Figure  8.  Data  points 
are  omitted  for  clarity.  Differences  a^e  shown  in  Figure  8  between  annealed 
and  unannealed  copper,  indicating  strain-hardening  effects  and  between  static 
and  impact  data,  indicating  strain-rate  effects. 

The  data  for  aluminum  in  Figure  2  show  some  effects  of  friction  or 
strain  hardening  and  shew  a  large  strain-rate  effect.  This  is  shown  by  the 
difference  in  position  of  the  data  points  from  static  and  impact  tests. 

In  Figures  3  to  7 ,  both  static  and  dynamic  data  concerning  five 
different  plastics  are  shown,  in  each  one,  a  strong  rate  dependency  is 
exhibited  showing  that  it  is  stronger  under  impact  conditions  than  under 
static  conditions.  It  is  also  shown  that  vinaar  impact  conditions,  an  apparent 
" s  rain- softening"  effect  exists  for  Nylon  and  polypropylene. 

A  quantitative  analysis  of  these  data  has  not  yet  been  made  to 
determine  the  coefficients  in  Equation  (14)  or  the  coefficient  of  friction  of 
Equation  (15). 


2 . 3  Dynamic  Measurements 

An  experimental  arrangement  for  measuring  the  motion  of  the  disk 
being  crushed  is  shown  in  Figure  9. 


lO 

t'v 


O 

lO 


LO 


o 

o 


LO 

t''. 


o 

ID 


LO 


<0 

Lh 

\ 

w 

0) 

«— I 

3 

o 


a 

w 

w 

(0 

s 

4-> 

f0 

,V 

cn 

Q 

*4-1 

o 

>4 

O' 

1_ 

a) 

c 

u 

u 


o 

a) 

a 

C/j 


fO 

4 — * 

"O 

d 

w 

a 

o 

4-J 

V-* 

<0  o 

Q.  *4—4 

o  w 

o  s 

L- 

3 

.  o 


o 

4—* 

f0 

1— 

fU 

a 

d 

w 


£ 

P 

o 


00 

d 

1-. 

3 

O' 


d 

d 

CO 


o  o 


*44 

N 


N 


-10- 


-11- 


Specific  Energy  of  Disk  Material,  E  (Joules/gram) 


1.0 

0.9 

0.8 


0  Dynamic 
A  Static 


-13- 


Fixed  Anvil 


Light 


Movable  Anvil 

- - >- 


Slit  of  light 


^  Disk 


Side  View  of  Anvils  and  Disk 


FIGURE  9.  Experimental  Arrangement  for  Measuring  Displacement 
and  Velocity  of  Moving  Anvil. 


An  oscillogram  of  the  light  intensity,  which  is  proportional  to  distance  of 
separation  of  anvils,  and  the  differential  of  this  signal,  which  is  proportional 
to  velocity,  vs,  is  shown  in  Figure  10. 


FIGURE  10.  Displacement  and  Velocity  of  the  Moving  Anvil 
of  Figure  9. 
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A  plot  of  the  acceleration,  obtained  by  graphical  differentiation  of  the  dis¬ 
placement  and  velocity  curves,  is  shown  in  Figure  11.  The  data  for  this  shot 
are  given  in  Table  I. 

The  energy  going  into  the  disk,  computed  from  the  dynamic  data, 
agrees  closely  with  that  obtained  from  measurements  of  rod  velocity  and  mass. 
The  low  rebound  velocity  of  the  rod  indicates  that  less  than  5  per  cent  of  the 
kinetic  energy  of  the  rod  goes  into  recoverable  elastic  energy  of  anvils  and 
disk  at  these  impact  velocities. 

Data  for  a  si  .lilar  shot  at  higher  velocity  are  shown  in  Table  II. 

A  rough  estimate  of  the  ratio  of  power  going  into  the  acceleration 
of  material  to  that  being  dissipated  in  the  material  is  given  by  comparing  the 
rate  of  change  of  kinetic  energy  of  the  material  with  the  rate  of  change  of 
internal  energy.  These  energies  are  plotted  in  Figure  12  for  the  data  of  Table  I. 
Note  the  different  scales  for  internal  energy  and  kinetic  energy.  The  maximum 
slopes  indicate  that  the  maximum  power  going  into  internal  energy  is  60  times 
that  going  into  kinetic  energy.  Impact  velocity  would  have  to  go  up  from  3.52 
x  10^  cm/sec  to  about  2.8  x  10^  cm/sec  for  power  going  *nto  kinetic  energy 
to  equal  that  going  into  internal.  Unfortunately,  tests  cannot  be  made  to  see 
how  internal  energy  behaves  at  this  impact  velocity  since  the  disk  would  be  so 
thin  that  friction  with  the  anvil  faces  would  dominate  the  process;  also  the 
anvils  could  not  be  considered  rigid. 


As  yet,  the  dynamic  measurements  have  been  used  only  to  compare 
integrated  power  data  with  total  impact  energy  computed  from  rod  velocity  and 
mass  and  to  compare  instantaneous  pressure  with  static  testing-machine 
pressure  data.  No  effort  has  been  made  to  find  the  time  function  which 
satisfies  the  equations  of  motion.  It  is  interesting  to  note  that  the  maximum 
pressures  indicated  by  the  two  dynamic  shots  are  7.9  x  10®  and  8.6  x  10® 
dyne/cm^,  'Hie  agreement  between  these  two  shots  at  different  velocities 
allows  a  good  estimate  to  be  made  of  the  strength  of  copper  under  flow  conditions. 
A  disk  pressed  to  similar  thickness  in  a  testing  machine  shows  higher  pressure, 
about  12  x  10®  dyne/cm^.  The  machine-tested  disks  were  barrel  shaped  after 
compression  while  those  tested  at  high  velocity  were  nearly  right  cylinders 
with  straight  sides.  This  probably  indicates  greater  friction  in  the  testing 


machine. 
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DATA  FROM  ROD- DISPLACEMENT  MEASUREMENTS  FOR  IMPACT 
OF  A  HARDENED-STEEL  ROD  ON  A  COPPER  DISK 
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FIGURE  12.  Kinetic  Energy  of  Disk  Material  Compared  with  Internal 
Energy  of  Material  for  the  Data  of  Table  I, 


3.  CONCLUSIONS 


The  disk-crushing  experiment  gives  quantitative  measurements  of 
the  energy  involved  in  the  flow  of  solid  materials  and  is  useful  where  the 
material  to  be  tested  is  less  rigid  than  the  anvils.  The  assumptions  made  in 
the  simplified  analysis  are  adequate  ter  small  deformations  and  can  be  corrected 
from  the  experimental  results.  Dynamic  measurements  can  be  made  and  the 
effects  of  strain  rate  and  st  ain  hardening  on  the  strength  of  materials  can  be 
determined. 


These  data  allow  ar>  immediate  means  of  comparing  materials  as  to 
their  resistance  to  flow  deformatior .  They  provide  necessary  data  for  any 
complete  theoretical  analysis  of  the  impact  problem. 


■  2A- 


4.  FUTURE  WORK 


The  disk-crushing  experimental  work  will  be  extended  somewhat. 

The  plastics,  which  have  strain-rate  dependent  strength,  will  be  tested  at 
intermediate  strain  rates  using  a  gun  or  a  drop  tower.  The  strength  of  frangible 
materials,  such  as  fiberglass-reinforced  plastics  or  acrilics,  will  be  tested 
using  the  dynamic  methods  described,  since  before-and- after  measurements 
are  not  possible.  The  experiments  on  friction  will  be  analyzed  and  extended. 
The  purpose  of  these  tests  is  to  make  quantitative  measurements  of  the  pnysical 
properties  of  materials  under  flow  conditions  and  to  express  these  properties 
by  appropriate  equations  such  as  Equation  (14). 

Shear  and  tension  properties  will  be  measured  by  shooting  through 
thin  sheets  or  into  wires  of  the  material  being  tested.  This  is  an  extension 
of  work  reported  previously.  The  relationships  between  deformation,  energy 
dissipation  in  the  material,  and  projectile  energy  loss  will  be  determined. 

The  compression,  shear  and  tension  data  will  be  used  to  formulate 
simplified  mathematical  model?  of  the  impact  problem  and  to  rate  pure  materials 
and  composites  as  to  their  method  of  failure,  the  energy  involved,  and  their 
ability  to  withstand  penetration. 

The  computer  program  reported  earlier  will  be  simplified  in  an  effort 
to  obtain  earlier  limited  results.  The  unique  predictor- corrector  scheme  used 
previously  will  be  developed  further.  The  tensor  analysis  of  deformation  will 
be  severely  simplified  and  modified.  The  easily  measured  material  properties 
coming  from  the  disk-crushing  experiment  and  from  the  tensile-failure  tests 
will  be  used  rather  than  a  tensor  formulation  of  material  properties. 
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